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o 

(N- We introduce a one-dimensional two-component system with the self-focusing cubic nonlinear- 

, ity concentrated at a symmetric set of two spots. Effects of the spontaneous symmetry breaking 

^ ■ (SSB) of localized modes were previously studied in the single-component version of this system. 

■ In this work, we study the evolution (in the configuration space of the system) and SSB scenarios 
for two-component modes of three generic types, as concerns the spatial symmetry of each com- 
ponent: symmetric-symmetric (Sm-Sm), antisymmetric-antisymmetric (AS- AS), and symmetric- 
antisymmetric (S-AS) ones. In the limit case of the nonlinear potential represented by two 5- 
functions, solutions are obtained in a semi-analytical form. They feature novel properties, in com- 

, parison with the previously studied single-component model. In particular, the SSB of antisymmet- 

■ ric modes is possible solely in the two-component system, and, obviously, S-AS states exist only in 
the two-component system too. In the general case of the symmetric pair of finite-width nonlinear 
potential wells, evolution scenarios are very complex. In this case, new results are reported, first, 
for the single-component model. These are pairs of broken-antisymmetry modes, and of twin-peak 
symmetric ones, which are generated by saddle-mode bifurcations separated from the transforma- 
tions previously studied in the the single-component setting. With regard to these findings, complex 
scenarios of the evolution of the two-component solution families are realized in terms of links con- 

^ ■ necting pairs of modes of three simplest types: (A) two-component ones with unbroken symmetries; 

QQ \ (B) single-component modes featuring density peaks in both potential wells; (C) single-component 

^SJ ■ modes which are trapped, essentially, in a single well. 

^\ 

^ . I. INTRODUCTION AND THE MODEL 

o 

A fundamental effect caused by the interplay of nonlinearity with symmetric potentials is spontaneous symmetry 
breaking (SSB). The simplest setting where this effect occurs is represented by double- well potentials. A commonly 
known property of one-dimensional quantum mechanics is that the ground state shares the symmetry of the underlying 
double- well potential If the self-attractive cubic nonlinearity is added to the consideration, the corresponding 
Schrodinger equation is transformed into the Gross-Pitaevskii equation for a Bose-Einstein condensate (BEG) loaded 

■ into the double- well potential [2j], or the nonlinear Schrodinger equation for photonic counterparts of the system 
0. The nonlinear term may break the symmetry of the ground state, replacing it by an asymmetric one which 
minimizes the energy of the system, provided that the strength of the nonlinearity exceeds a certain critical value (see 
Ref. (4,] for the general analysis, and Ref. [5,] for the consideration of the SSB of self-trapped states in BEG). The 
experimental realization of the SSB in double- well potentials was reported in BEG @ (this was done in the condensate 
with the self-repulsive nonlinearity, which implies the spontaneous breaking of the antisymmetry of the lowest-energy 
antisymmetric mode) and in nonlinear optics [7|. 

The symmetry-breaking bifurcation, which destabilizes the original symmetric ground state and gives rise to the 
SSB in the nonlinear systems, was originally predicted in a discrete model of self-trapping 0]. In nonlinear optics, a 
similar bifurcation was analyzed in Ref. @ for continuous-wave (spatially uniform) light si gna ls in dual-core fibers. 
The allied bifurcation for solitons in nonlinear dual-core fibers was studied in detail in Ref. [lOj. Later, the SSB was 
studied for gap solitons in the model of the dual-core fiber Bragg gratings with the same cubic nonlinearity as in the 
ordinary fibers [ll|. The SSB effects were also predicted for matter- wave solitons in the self-attractive BEG loaded 
into a double-channel potential trap [l^-fTsj. 

Typically, the self-attractive cubic nonlinearity gives rise to the soliton bifurcations of the subcritical (backward) 
type (in other words, these are phase transitions of the first kind). In that case, branches of asymmetric modes 
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emerge as unstable ones, going backward (i.e., in the direction of the decrease of the soliton's norm (energy)). They 
get stabihzed after switching the evolution direction forward at turning points of the bifurcation diagram [16| . It was 
also demonstrated that the combination of the self-focusing nonlinearity with a periodic potential applied in the free 
direction (perpendicular to the direction of the action of the double-well potential) tends to change the character of 
the soliton bifurcation from sub- to supercritical (i.e., replace the phase transition of the first kind by a transition 
of the second kind). In the supercritical case, the asymmetric branches emerge as stable ones, going in the forward 
direction, i.e., in the direction of the increase of the soliton's norm [l^, The bifurcation for gap solitons in the 
dual-core fiber Bragg grating is of the forward type too [llj . 

A challenging alternative to the use of the ordinary linear double-well potential is the setting based on an effective 
nonlinear potential, alias pseudopotential (as it is called in solid-state physics [l7|). which is induced by the double-peak 
spatial modulation of the local strength of the self-focusing nonlinearity. Pseudopotentials featuring the double-well 
shape can be realized in optics and matter waves [isl [l9| . The ultimate form of such a setting is the one with the 
nonlinearity concentrated at two points, represented by a symmetric pair of delta-functions or narrow Gaussians 
[20I I21I I (the simplest model of that type, with the self-attractive nonlinearity represented by a single delta-function, 
was first introduced in Ref. p^l. Two-dimensional counterparts of the system were considered too, in the form of 
two parallel channels [1^, or two symmetric circles [1^, in which the self- attractive nonlinearity is localized . The 
SSB of solitons in the symmetric nonlinear double- well potentials was studied in detail, featuring the subcritical type 
of the symmetry breaking H^, HH, [111 . 

It is relevant to mention that SSB effects were also analyzed for discrete solitons in dual-core nonlinear chains, with 
the uniform linear coupling between them (27j , as well as with the linear coupling which links a single pair of sites in 
the parallel chains [l^ . The form of the linear coupling determines the type of the respective SSB bifurcation, which 
is subcritical in the former case, and supercritical in the latter situation. The action of the nonlinearity concentrated 
at a pair of sites embedded into a single linear chain was studied too [2^j30|. The SSB effects in this setting (for both 
straight and circular host linear chains) were analyzed in recent work [3o|T^ The SSB in a similar system, with a pair 
of symmetric nonlinear sites side-coupled to the infinite linear chain, was recently studied in Ref. [31| . 




An obviously relevant generalization of the work outlined above is extension to two-component systems. In BEG, 
this implies a mixture of two different hyperfine states of the same atomic species [2] , while in optics it implies the co- 
propagation of light beams mixing two different polarizations of light, or different carrier wavelengths. Recently, the 
SSB, along with the related dynamical effect of Josephson oscillations between the wave functions trapped in adjacent 
potential wells, have been analyzed in diverse models of binary systems embedded into ordinary (linear) potentials 
of the double- well type [S^. An experiment was performed too, for the corresponding BEG mixture [S^. However, 
SSB effects have not yet been studied in two-component systems trapped in double- well nonlinear (pseudo)potentials. 
A straightforward possibility to implement the latter setting can be found in nonlinear optics, where the spatially 
nonuniform nonlinearity, which forms the double- well structure, will act similarly on both polarization components, in 
the form of the SPM (self-phase modulation), as well as on the XPM (cross-phase-modulation) nonlinear interaction 
between them @ . 

In this work, we aim to perform the analysis of the two-component system in the nonlinear double-well potential, 
with emphasis on manifestations of the SSB in the two-component mixture. To this end, we aim to consider the basic 
setting with the symmetric pair of strongly localized nonlinear spots embedded into the one-dimensional linear host 
medium. Thus, our model is the generalization, for wave functions <f){x,t) and ^(x,t) of the two components, of the 
single-component model introduced in Ref. 20] : 



where G > is the relative strength of the XPM nonlinearity, while the SPM coefficient is normalized to be 1. In the 
optical model corresponding to these equations, the evolutional variable t is actually the propagation distance (usually 
denoted as z), while x is the transverse coordinate (in the BEG model, t is simply time). As for the XPM coefficient, 
its typical values in optics are G — 2/3 for the coupled linear polarizations, and G = 2 for the pair of orthogonal 
polarizations or different carrier wavelengths [1]. Other values of G are also possible, for elliptically polarized beams. 

Following Ref. f20i], the nonlinearity- modulation function in Eqs. ([IJ and ([2]), which corresponds to the symmetric 
set of two strongly localized nonlinear spots, is adopted in the form of 





(1) 



(2) 




(3) 
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which is subject to the normahzation condition, 



+00 



g(x)dx 



(4) 



Profiles of modulation function ([3]), which keeps the double- well structure at a < \/2, are shown, for different values 
of a, in Fig. [TJ 




Figure 1: (Color online) Shapes of the double-well nonlinearity-modulation function ((3}, normalized to its maximum value, are 
shown (as per Ref. ^2d\) for different values of scaled width a of the individual well. 



Equations ([T]) and ^ conserve three dynamical invariants, viz., two norms and the energy (Hamiltonian), 



M 



H = 



+ 00 



\(i)(x)\^dx, N 



+ 00 



m^)?dx, 



+ I^.l' + (101' + IVI' + 2G 101' IVf 



dx. 



(5) 
(6) 



Stationary localized solutions to Eqs. ([T]) and ^ are sought for as 4) — e~'^^*u{x), ip = e~*^*w(a;), where the 
chemical potentials of localized modes must be negative, A < 0, < 0, and real functions u{x) and v{x) satisfy 
equations 



Xu + -u" H ■j= 

2 ClJiT 



tJ.v + -v' 



1 



exp 



exp 



(x + iy 



(x + iy 



exp 



exp 



{x-lf 



(x-iY 



{v^ + Gu^)v = 



A special version of the model corresponds to the limit of a -^• 0, with Eqs. ([7]) and ([5]) going over into 

Xu = -{l/2)u" ~[6{x + l) + 6{x~l)]{u^ + Gv^)u, 



fJiV 



-{l/2)v" - [5{x + 1) + 5{x - 1)] + Gu') V, 



(7) 
(8) 



(9) 
(10) 



where 6{x) is the Dirac's delta-function. 

To conclude the Introduction, it is relevant to mention that a two-field model with a symmetric set of two strongly 
localized nonlinear spots was recently introduced for the quadratic (second-harmonic-generating) optical nonlinearity 
SJ). In the limit case when the nonlinearity profile is described by the pair of (5- functions, cf. Eqs. ([9|) and (|10p . 
the corresponding equations for the complex amplitudes of the fundamental-frequency and second-harmonic waves, 
U{x, z) and V {x, z), take the form of 
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iU, + {l/2)U^, + S{x)U*V = 0, 



(11) 



2iV, + il/2)V^,-QV + il/2)6ix)U^ = 0, 



where the asterisk stands for the complex conjugate, z is the propagation distance, and real constant Q measures the 
mismatch frustrating the parametric interaction between the two waves. Although involving the two wave fields, this 
system ffor which exact solutions for symmetric, antisymmetric, and asymmetric localized modes have been found in 
Ref. |34| ) is actually closer to the single-component model in the case of the cubic nonlinearity, due to the coherent 
character of the nonlinear coupling between the waves U and V. In particular, system (jlip directly reduces to its 
single-component cubic counterpart in the cascading limit, which corresponds to large positive values of Q. 

The rest of the paper is organized as follows. The limit case of the two-component system with the modulation 
profile based on the set of (5- functions, which corresponds to Eqs. © and (|10p . is considered in detail in Section 2 (in 
the beginning of the section, we also give an exact solution for the two-component system with the single nonlinear 
(5- function, as this simple case was not considered previously) . The system with the pair of ^-functions admits a semi- 
analytical solution, which reduces to a system of four coupled cubic algebraic equations for amplitudes of asymmetric 
modes, while symmetric (Sm) and antisymmetric (AS) ones can be found in a fully analytical form. 

In Section 3, we report results of the analysis of the two-component model in the general case, corresponding 
to the nonlinearity-modulation function ([3]) with finite a. In that case, the analysis is based on a comprehensive 
numerical solution of ODEs ([7]) and © with a > 0. The consideration in Section 3 starts by revisiting the single- 
component model, with the aim to report additional solutions of the Sm and broken-AS types (the latter means 
modes with broken antisymmetry). These solutions were not reported in Ref. [20|, as they represent pairs of modes 
generated by saddle-node bifurcations which occur at sufficiently large a, being detached from the branches which 
continuously evolve from the limit case of a — > 0. The subsequent analysis of the two-component system reported 
in Section 3 produces a full picture of bifurcations for compound modes of the Sm-Sm (symmetric-symmetric), AS- 
AS (antisymmetric-antisymmetric), and S-AS (symmetric-antisymmetric) types, and their broken counterparts, i.e., 
modes with spontaneously broken (anti)symmetries. Because the full picture turns out to be extremely complex, we 
report only parts of it, which may be presented in a relatively simple form. Despite the complexity of the picture, 
some underlying principles can be formulated in a general form: the fundamental role is played by simplest solutions, 
namely the two-component unbroken modes (i.e., those with the unbroken (anti)symmetries), and single- component 
solutions of both the unbroken and broken types. As a matter of fact, numerous families of two-component solutions, 
which feature sophisticated evolution in the system's configuration space, may be understood as branches linking 
various pairs of these simplest modes (the complexity is often produced by the fact that the branches make one or 
several loops, in the course of their evolution). 

The paper is concluded by Section 4. In particular, we discuss the necessity of the systematic analysis of the 
dynamical stability of the diverse stationary modes found in this work. With a few exceptions, we do not analyze the 
stability here, as the study of the stationary modes and their bifurcations is by itself a sufficiently heavy topic for a 
single paper. 



Before presenting results for the model with the symmetric pair of (5-functions, it makes sense to briefly consider 
the case of the nonlinearity represented by the single (5-function, when Eqs. (H]) and © are replaced by 



II. THE MODEL WITH THE DELTA-FUNCTIONS (a ^ 0) 



A. The single delta-function 




(12) 



(13) 



and, accordingly, their stationary version, Eqs 



^ and (Uni), is replaced by 




(14) 
(15) 
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It is straightforward to find solutions to Eqs. p4)) and (115]) with both fields different from zero: 

{u{x),v{x)} = {A, B} exp(~^/~2{X,fi}\x\] 
V2 



{A^B^} = ^l—[G^^{fIJ^-^^{X;fI} 



(16) 
(17) 



These solutions exist for G 7^ 1 (G = 1 is the degenerate case corresponding to the Manakov's nonlinearity Q), for 
negative chemical potentials A and /i, in the regions of 



(18) 
(19) 




G < W - < if G< 1. 
fi G 



It is also instructive to cast Eq. p7|) into the form of relations between the chemical potentials and norms M and 
N of the two components (see Eq. ©). Obviously, {M,N} = {A^,B^} /^-2{A,^} for solutions (HH), hence 



{M,N} = 




(20) 



Although we chiefly leave the issue of the dynamical stability aside in this paper, Eq. (I^Hl) makes it possible to 
address the stability problem, using the Vakhitov-Kolokolov (VK) criterion. It states that, for the single-component 
model, with the soliton family characterized by dependence AI = M(A), the necessary stability condition is dAI/dfj, < 
36] , while for the two-component model with separately conserved norms and two independent chemical potentials 
the switch between the stability and instability occurs at points where the corresponding Jacobian determinant 
vanishes: 



d (M, N) 



a(A,M) 



0. 



(21) 



In the single-component model with the single nonlinear (5-function, the soliton family is degenerate^ in the sense 
that its norm does not depend on the chemical potential: M = 1, hence this family seems as a neutrally-stable 
one, which actually does not guarantee the true stability. In reality, the single-component family was found to be 
completely unstable [21]. In the present case, the calculation of Jacobian J {X, fj.) for the two-component solitons per 
Eq. (I20I) demonstrates that this determinant identically vanishes, J (A, /i) = 0, suggesting that the total family of 
solutions (|16p. ()17p is completely unstable too. Indeed, one can easily find the following set of exact nonstationary 
solutions to Eqs. (fT2l) and (fT3l) . with arbitrary real constant a and /?: 



{(I){x,t)ij{x,t)} = {a,/3}|t|-i/2exp 

{a,b} = -{{a\l3^}+G{fi\a'})sgn{t), 



^^{\x\-i{a,b}f 



(22) 
(23) 



which generalize similar exact single-component solutions found in Ref. [21|. At i < 0, solutions ([22|) . (p3)) describe 
the collapse of the perturbed soliton (the formation of a singularity at x = 0), which happens at i — > —0, while at t > 
the solution describes decay of the soliton. Thus, the nonstationary solutions demonstrate that the solitons belonging 
to the degenerate family ([T6|. (ITTl) . supported by the single 5- function, suffer either the collapse or decay, being 
indeed unstable, as conjectured above. On the other hand, we expect that the replacement of the ideal (5-functions by 
a finite-width profile may stabilize the solutions, cf. Ref. [s^ where this was demonstrated for the single-component 
model. 



B. Amplitude equations for the system with two delta- functions 

Proceeding to the model with the symmetric set of two (5-functions, based on Eqs. © and ((TU]), our first objective is 
to obtain analytical results, following the pattern of Ref. [13], where a full analytical solution was found for stationary 
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states in the corresponding single-component model. Off points x = ±1, equations ([9]) and (|T0)) are linear, a general 
solution to which, decaying at |a;| — >■ oo, can be sought for as 



BieV2W(=^+i), at X < -1, 
u{x) = I Aoe-V^\(^-^) + BoeV2|^(^+i), at - 1< .x < +1, 
>lie-V2|A|(^-i), at X > +1, 



(24) 



v(x) — < 



DieV^\(^+^\ at X < -1, 
Coe-N/2H(^-i) + DoeV2H(^+i)^ at - 1< x < +1, 
Cie-V2H(^-i)^ at X > +1, 



(25) 



with constant amplitudes Aq, Ai, Bq, Bi and Co, Ci, Dq, Di. The continuity of the wave functions at x = ±1 imposes 
four relations on the amplitudes. 



Bi = Bo + Aoe^V^, A^ = Ao + Soe^V^, 
Di ^ Do + Coe^^^, Ci = Co + Doe^V^^. 



(26) 
(27) 



Using these relations, one may eliminate "inner" amplitudes Ao,-Boi Cq, i^o in favor of their "outer" counterparts 
^1 , Ci , £>! : 



^0 



Bq = 



=V2|A| _ 



<-^o — . /—— , ^0 — 



.W^M - 1 



eV2|A| _ 1 
e^v^Ci - £>! 



(28) 
(29) 



Further, the integration of Eqs. (jO]) and (jlOp in infinitesimal vicinities of points x = ±1 yields expressions for jumps 
(A) of the first derivative at these points, 



A(u') U=±i 
A(i'')U=±i 



(30a) 
(30b) 



The substitution of solution ([M)) and psp into relations (pO| leads to a system of equations for the outer amplitudes, 

^Bi -Bo + Aoe^V^^ 



= {Bf + GDl)Bi, 

^/\X\/2(^Ai-Ao + Boe^V^\) = {Aj + GC^,)A,. 

/H72 - + Coe^V^) = {Dl + GBf)Di, 

v4^(Ci-Co + i?oe2V^) = {Gf + GAl)Gi. 



(31) 
(32) 

(33) 
(34) 

After the substitution of expressions ([28]) and (f29l) into Eqs. ([32|) -(|34 |) . we end up with a system of four coupled cubic 
equations for Ai,Bi,Ci, Di: 



e^V^^Bi - Ai) = sinh (2^/2\X\^ {Bf + GDI) Bi, 
e^v^Ai - Bi) = sinh (2V2|A[) {Aj + GCf) Ai. 
g2y2H£,^ - Ci) = sinh (2v/2R) (ZJ^ + GBf) Di, 



(35) 
(36) 
(37) 



7 



[e^^^Ci ~ = sinh (2v/2R) {Cf + GAl) Ci. 



(38) 



Symmetric solutions to Eqs. ([35l) -(|38 l) (as said above, they will be referred to as "Sm-Sm" modes, i.e., those which 
are spatially symmetric in each component) are looked for by setting Ai = Bi and Ci — Di, which reduces the system 
to equations 



Al + GCf = 



C\ + GA\ 



l + e-V2|>l 



/2R 



l + e-V2IH 



(39a) 
(39b) 



Similarly, solutions to Eqs. (|35 |) - (|38|) which are antisymmetric in each component ("AS- AS" states) are looked for 
with Ax = -Bi and Ci = -Di, reducing Eqs. dM])-® to 



Al+GCf = 



Ci + GAi 



'2\X\ 



1 _ e-V2|A| ' 

1 _ g-2V2H ' 



(40a) 
(40b) 



There also exist solutions to Eqs. (|35 p -p8 p with mixed symmetry, e.g., antisymmetric in component u and symmet- 
ric in u, which implies Ai = — i?i and Ci — Di ("S-AS" states). This restriction reduces the system of Eqs. (|35l) - ([55)) 
to a simplified one, 



Al + GCf = 
Cl + GA\ = 



1 _ e-2V2|A| ' 



(41a) 
(41b) 



C. Analysis of the solutions 

Systems (|39|) . (|40p and (|4T|) are obviously analytically solvable, as each one is tantamount to a system of two 
inhomogeneous linear equations for Af and C^, while the general algebraic system p5l) - ([55|) can be only solved in 
a numerical form. Norms ([5]) and energy ([6]) for the solutions were also calculated numerically, because analytical 
expressions for these integrals are very cumbersome (they are actually analytically intractable even in the single- 
component model |20j). The results are displayed below for G — 2. As said above, in optics models this value 
corresponds to the interaction between orthogonal circular polarizations, or between waves with different carrier 
wavelengths. 

The numerical solutions reveals a plethora of different branches of the localized modes, with different degrees of 
the asymmetry. Typical examples of the corresponding spatial models of the different types are presented in Fig. [5] 
In this paper, we aim to report only a part of the results, which represent the most characteristic findings, that may 
be presented in a sufficiently clear form, as the full description of the (anti)symmetry breaking in the two-component 
model (both for a — > and finite a) is extremely complex. Generally, the results may be understood as those pertaining 
to modes generated by the spontaneous (anti)symmetry breaking from the above-mentioned solutions of three special 
types, i.e., Sm-Sm, AS-AS, and S-AS ones. This is the principal difference from the single-component model, where 
asymmetric modes could be generated solely from the symmetric ones (antisymmetric states were dynamically unstable 
at small a and stable at larger a, but they never underwent a symmetry-breaking bifurcation). 

Each species of the modes may be characterized by their energy ^ and spatial asymmetries of the two fields, which 
are defined as follows: 

(+00 \ 

I dx V(4)l' - j dx mx), V.(x)}|2 . (42) 
-oo / 

Exacerbating the complexity of the results, there is, in general, no one-to-one relation between chemical potentials 
and norms M,N, as well as total energy H, see Eqs. (O, (jH]). However, the asymmetries and energy of the solutions 
smoothly depend on M and N, allowing us to classify the solutions according to the type of their (a)symmetry 
(Sm-Sm, AS-AS, S-AS, and their "broken" versions). 
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Q> f * * 




Figure 2: (Color online) Typical examples of the exact modes in the model with two 5-functions, produced by solution of 
algebraic equations (|35p - (|38p at A = 1.5 and jU = 1.0: Sm-Sm modes (a), AS-AS modes (b), and S-AS modes (c), each with 
broken and unbroken (anti)symmetry. 



1. Solutions obtained from the symmetric-symmetric modes (Sm-Sm type) 

We start by the presentation of the results for modes generated by the SSB of the states of the Sm-Sm type, as this 
case admits direct comparison with the results obtained for the SSB of the symmetric mode in the single-component 
model (20l |. Examples of the symmetric modes with the broken and unbroken symmetry are shown in Fig. EJa), for 
A = 1.5 and /i — 1.0. In Fig. [3l the asymmetry of one component of the "broken" solutions of the Sm-Sm type is 
displayed, along with total energy ([B]), considered as functions of norms ([S]). The asymmetry of the second component 
actually follows the same pattern. 

A characteristic feature of the solutions of the Sm-Sm type, both "broken" and "unbroken" ones, is that they allow 
to establish a one-to-one correspondence between the norms and chemical potentials of the two components. It can 
be checked that both the norms and total energy ([6]) of these modes grow monotonously with the increase of |A| and 
\fj,\. It is worthy to note too that the energy of the symmetry-broken Sm-Sm modes (which is shown in Fig. El^b)) 
is always higher than the energy of their unbroken counterparts, for the same values of M and TV, which suggests 
that the "unbroken" and "broken" states may be stable and unstable, respectively, similar to the situation in the 
single-component model (20| . 

The cross section of the plots in Fig. [3] along M = and = correspond to the single-component model. 
This fact explains why the symmetry-broken modes exist only in a limited domain in the {M, N) plane, as in the 
single-component model with a — > the solutions with the broken symmetry exist in a narrow interval of values of 
the norm [20|], 

1< Af < 2/3 + (8/27) (3/4 + In 2) w 1.09. (43) 

The most salient feature observed in Fig. ^a.) is the lacuna ("hole") in the existence plot of the modes. In fact, the 
cross section of the plot, at a fixed value of M or N, which does not cut the lacuna (including the sections running 
through M — and = 0), represents a situation qualitatively similar to that observed in the single-component 
model: the asymmetry parameter (|42p varies continuously from —1 to -1-1, passing the zero value, while the norm takes 
values in a narrow interval, cf. Eq. P5)) . On the other hand, the cross section which cuts through the lacuna reveals a 
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qualitatively new situation, that does not occur in the single-component model: the asymmetry parameter varies in two 
disjoint outer intervals (which are, naturally, mirror images of each other), (AO^) . < |A00| < 1 (the corresponding 
interval of the variation of the norm remains narrow), while the inner segment, — (AO^) . < A00 < + (AO^) . , 
remains empty. Because there is no solution with AQ^ = in the latter case, this means that the corresponding 
bifurcation in the single-component model (if such a bifurcation were possible) would be completely different from 
the one observed in the actual single-component equation: it would be a bifurcation loop disconnected from the line 
of the solutions with the unbroken symmetry. 

Nevertheless, the analysis of the full results for the two-component system based on the set of two (5-functions 
demonstrates that branches of all the species of broken-symmetry solutions (those of the Sm-Sm, AS- AS, and S-AS 
types, see results for the two latter species below), which are parameterized by either chemical potential, A or /i, while 
the other one is fixed, actually originate, at bifurcation points, from the families of allied solutions with unbroken 
(anti)symmetry, and eventually merge back into the same families (see more details below). In other words, the 
bifurcation curves are indeed shaped as loops, but they are not detached from the respective (anti)symmetric solution 
families, which actually constitute parts of the bifurcation loops. 



2. Solutions obtained from the antisymmetric- antisymmetric modes (AS-AS type) 

As mentioned above, the single-component model with a — >■ does not give rise to spontaneous antisymmetry 
breaking of antisymmetric modes [l^l • In the present two-component system, double-antisymmetric modes (the AS- 
AS type) do admit SSB, which is therefore a novel effect, in comparison with the single-component model. In Fig. |4l 
the global picture of the AS-AS solution family with the broken antisymmetry is displayed by means of two asymmetry 
measures (H^ . AO^ and A0^, along with total energy ©, shown as functions of norms M and A^. In fact, the plots 
presented in panels (a) and (b) are tantamount to each other, but, being shown with respect to the fixed frame of 
the axes {M and N ), they provide mutually complementary views. To clarify qualitative aspects of the picture. 
Fig. [5] additionally displays the projection of the plots in Fig. |ll[a),(c) onto the planes of (M, AG^) and (M, IniJ), 
respectively. Comparison of Fig. \M,c) with Fig. Elb) clearly demonstrates that the AS-AS solutions have much higher 
energy than their counterparts of the Sm-Sm type, hence the AS-AS modes may be more vulnerable to instabilities 
(as their AS counterpart in the single-component model (20j). 

In connection to the pictures displayed in Figs. |3]and[Sl it is relevant to mention that, in the single-component 
model, the mode with the unbroken antisymmetry exists in an infinite domain limited from below, 1.84 < M < oo. 
On the contrary, we here observe that the norms of the "broken" AS-AS solutions take values in a limited domain. As 
seen in Fig. [5j the corresponding existence plot also contains a lacuna, and, in addition to that, it is confined to values 
jAO^I < 1. A formally similar bifurcation in the single-component model, should it be possible, might again seem as 
a loop detached from the unbroken-antisymmetry modes, which does not reach the limit values of jAG^^^j = 1. 

To consider the bifurcation which actually accounts for emergence of the "broken" AS-AS solutions from their 
counterparts with the exact antisymmetry, it is convenient (as said above) to fix one of the chemical potentials, A 
or fi, varying the other one. The analysis of the numerical results demonstrates that the antisymmetry-breaking 
bifurcation occurs, with the variation of A, when takes values exceeding (|/^|)^i„ = 0.0601 (indeed, for — >• 
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(a) (b) 




Figure 4: (Color online) AS-AS modes with broken antisymmetry. Panels (a,b) and (c) display, severally, the asymmetries of 
components and ip, defined as per Eq. ([42}, and the total energy vs. norms M and A'' of the two components. 




Figure 5: (Color online) The projections of the three-dimensional plots from Fig. |4ja) (left panel) and Fig. [U^c) (right panel) 
onto the two-dimensional planes of (M, AO^) and {M,\nH), respectively. 

one should actually return to the single-component model parameterized by A, where the AS mode does not undergo 
any bifurcation). Emerging from the branch with the unbroken antisymmetry, the solution of the broken AS-AS 
type eventually merges back into the original branch, as |A| attains a certain maximal value (which increases with 
the increase of fixed |/i|). If |A| and/or remain small enough, the numerical results demonstrate a one-to-one 
correspondence between the chemical potentials and the corresponding norms, which grow monotonously with the 
absolute values of the chemical potentials. However, unlike the modes of the Sm-Sm type, for the AS-AS solutions 
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this correspondence breaks down at larger values of |A| and 

If the fixed chemical potential exceeds another critical value, — 0.198, an additional bifurcation occurs on the 

broken- AS- AS branch, giving rise to two additional solution families. With the further increase of the absolute value 
of the varying chemical potential, |A|, at some point |A| — {\M)rrLa.x '-'^'^ '^^ secondary branches merges with the 
original one from which it had emerged. The other secondary branch extends to larger values of |A|, and eventually 
merges back into the underlying family of the solutions with the unbroken antisymmetry. The range where the 
secondary branches are found expands with the increase of fixed \ fi\. Another generic feature of the bifurcation, which 
is responsible for the merger of the broken- AS- AS branch back into its unbroken counterpart, is that this happens 
when one norm is much larger than the other - typically, at M ~ 0.05, while iV ~ 2. 

3. Solutions obtained from the symmetric-antisymmetric modes (S-AS type) 

The existence of modes with unbroken or broken mixed symmetry, of the S-AS type, is another obvious difference 
of the two-component system from its single-component predecessor. As might be expected, these solutions exist 
only in the region of the (M, N) plane where the Sm-Sm and AS- AS states coexist. The global structure of the 
S-AS family with the broken symmetry and antisymmetry is displayed in Fig. [6] by means of plots showing the 
asymmetry measures of its initially symmetric ((f)) and antisymmetric (tp) components, as well as the total energy, 
vs. the norms. Since the pattern of the solutions is rather complex, we complement Fig. [5] by Fig. [71 where the 
projection of the three-dimensional plots of the asymmetry measure of components cf) and ip from Figs. [nia),(b) onto 
the two-dimensional planes of (M, AO^) and (M, A0^), respectively, is presented. 



(a) (b) 




Figure 6; (Color online) The family of mixed-symmetry S-AS modes, with the broken (anti)symmetry. (a) and (b): The 
asymmetry measures of the symmetric and antisymmetric fields, (p and tp, respectively, (c) The energy of the S-AS modes with 
the broken and unbroken (anti)symmetries vs. norms M and A'^. 

These plots are very different from those for the Sm-Sm and AS- AS modes, cf. Figs. |3][S] To consider the relation 
between the mixed-symmetry modes with the their Sm-Sm and AS-AS counterparts, it is relevant to stress, that, 
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Figure 7: (Color online) Projections of the three-dimensional plots from Fig. (left panel) and Fig. [5] (b) (right panel) onto 
the two-dimensional planes of (A/, AG^) and (A/, A0^), respectively. 

as may be naturally expected, the branch of S-AS solutions with unbroken (anti)symnietries smoothly merges into 
the corresponding solutions on the Sm-Sm type in the limit of the vanishing norm N in the AS component. On 
the other hand, when the latter norm (N) increases, the energy of the system increases too, and the configuration 
smoothly approaches the opposite limit of the AS- AS solution. If chemical potential |A| is large enough, the latter 
limit corresponds to vanishing norm AI of the Sm component. However, for |A| < 0.71, the A^f-norm does not vanish 
in the limiting case; instead, it approaches some finite value, whereas the A^-norm rapidly diverges, as seen in Figs, 
inija) and (b). In fact, this singular behavior is a peculiarity of the model with the (5-functions (a — > 0), while at a > 
(in the model with a finite width of the pseudopotential wells) the situation is different, as shown in the next section. 

As well as in the cases considered above, the states of the S-AS type with the broken (anti)symmetries emerge from 
their "unbroken" counterparts through a bifurcation (i.e., they are connected). We discuss this bifurcations in more 
detail below, considering the model with finite a. 

It can also be checked, using detailed results of the numerical analysis, that, similar to the case of the Sm-Sm modes, 
the norms and energy of the S-AS states monotonously grow with the increase of |A| and |/Lt|, which, in particular, 
means that the norms are in one-to-one correspondence with the chemical potentials, also similar to the case of the 
Sm-Sm modes, and in contrast to their AS- AS counterparts. Lastly, comparing Fig. [HJc) to Fig. |3|b), we conclude 
that, for the same values of norms M and N, the energy of the S-AS modes is essentially higher than the energy of the 
fundamental Sm-Sm states, which suggests that the S-AS states may be more vulnerable to dynamical instabilities. 

III. THE REGULAR MODEL (a > 0) 
A. Preliminaries 

The results reported in the previous section were based on the numerical solution of coupled cubic equations ([55)1 - 
(I38p , which were derived from the analysis of the singular version of the model, with the nonlinear potential represented 
by the symmetric set of two (5-functions, see Eqs. (0) and (ITUl) . It is known from the analysis of the single-component 
model that this limit case, corresponding to a — in Eq. ([3]), is degenerate in several aspects (in particular, the 
symmetric mode exists only up to a finite value of the norm, see Eq. (|43p ). In the single-component model, the 
degeneracy was lifted by the transition to the nonlinearity-modulation function in the general form (j3|), with finite 
a (recall that the modulation function keeps its double- well shape up to a = V2). Another generic difference of 
the situation with finite a is that, unlike the case of the nonlinear potential represented by the (5-functions, smooth 
modulation function ([3]) may support localized patterns with local maxima shifted from the bottom points, x = ±1. 

We have expanded the analysis of the two-component system to the general case of g{x) taken in the form of ([3]) 
with a > 0. As demonstrated below, the shape of various modes and bifurcations of the respective branches strongly 
alter, against the limit of a — > 0, starting from values a > 0.1 (in the single-component model, the essential change in 
the picture of stationary modes was observed in, roughly, the same range (20|). All the types of the two-component 
modes, i.e., Sm-Sm, AS- AS, and S-AS, in their symmetry-broken and unbroken forms alike, feature strong changes 
against the case of a —> in this region. In fact, these changes are strongly affected by additional modes with broken 
(anti)symmetry, which appear, at a large enough, still in the single- component model. These modes were not found 
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in Ref. [20j, because they are generated "from nothing" by saddle- node bifurcations, that are disconnected from 
the symmetry-breaking bifurcation for the symmetric mode, and the stabihzation transition for the non-bifurcating 
antisymmetric one, on which the analysis was focused in Ref. [13] ■ Therefore, in this section we first present the 
new results for the single-component model, and then report, in a rather brief form, changes which occur to the 
compound modes in the two-component system. A full description of the bifurcations in the two-component system 
with a G (O, turns out to be extremely complex, as there is a large number of ramifications of the bifurcation 
scenarios. 

Some comments are relevant too as concerns the numerical methods used for the search of localized solutions to 
ODE system ([7]), ([8]) with a > 0. In addition to the obvious fact that the numerical analysis had to cover the three- 
dimensional parameter space, (a, A, fj,), one of difficulties in the numerical study is related to the fact that, for some 
solutions that had to be obtained, components u{x) and v{x) feature very different (up to a few orders of magnitude) 
variation rates. Another problem is that, because for small a the coupling between the equations is confined to small 
neighborhoods of points x = ±1, it was necessary to use a numerical scheme with a higher density of grid points in 
those areas. 

One numerical method that was used to solve the system of equations ([7]) and ([5]) was based on the shooting 
algorithm [381] . which is, actually, better suited for producing the solution at smaller values of a. Another approach, 
using the centered-difference-approximation method, turns out to be more appropriate for larger a. In particular, 
the application of the shooting method was assisted by the fact that, for smallest values of a, the solution has to 
be close to that available in the analytical form given by Eqs. ([M)) and (P5|) for a — >■ 0. This fact provided a good 
initial guess for the shooting method at a < 0.4. An important issue was also the choice of the shooting point. The 
method did not work with these points taken at ±" infinity" , but it worked well when the point was set in the middle 
of either nonlinear-potential well (see Eq. ([3])). By means of this technique, 5 x 10^ solutions for the localized modes 
had been collected, which made it possible to formulate general conclusions about the bifurcation scenarios in the 
two-component system. 

B. New results for the single-component model 

In Ref. f23|, no bifurcations were found on the branch of single-component antisymmetric modes, both for a — >■ 
(where this result was exact) and for finite a (in the numerical form). Our analysis complies with those conclusions. 
Nevertheless, as shown in Fig. |51 our numerical studies, performed at finite a, reveal saddle-node bifurcations (more 
than a single one, see Fig. [5Uc,d)), that generate pairs of single-component modes with broken antisymmetry, which 
are disjoint from the branch with the unbroken antisymmetry (for this reason, they were not found in Ref. 120]). An 
example of the pair of the newly found solutions with the broken antisymmetry is displayed in Fig. ^sl). An essential 
peculiarity of these modes, which is seen, for instance, in the case of the lower mode displayed in Fig. ^a), is that 
they tend to be trapped in a single nonlinear-potential well, instead of being spread between both, which is the case 
for the mode with the unbroken antisymmetry. It will be shown below that such modes, essentially trapped in a single 
pseudopotential well, play an essential role in bifurcation scenarios of compound modes of the AS-AS and S-AS types 
in the two-component system at finite a. 

Furthermore, we observe similar saddle-node bifurcations of symmetric and asymmetric solutions in the single- 
component model. In particular, a pair of new symmetric-solution branches, generated by such a bifurcation, is 
displayed in Fig. [TU] These bifurcations take place on the "old" symmetric and asymmetric branches at some critical 
values of a (in contrast to the bifurcations displayed in Fig. [51 which are disconnected from the "old" branch of 
the antisymmetric solutions). For the asymmetric and symmetric solutions, this happens at a = 0.67 and a — 0.72, 
respectively (the bifurcation which gives rise to the pair of the asymmetric branches is not displayed here). As a grows 
beyond these critical values, the new branches move away from the parent ones. It may happen that more pairs of 
branches emerge through additional saddle-node bifurcations (cf. Fig. [51 where two such bifurcations are shown for 
the broken antisymmetric solutions), but we did not aim to look for them. 

The difference of the newly emerging symmetric states from the original one, which evolves continuously starting 
from a = 0, is that, at relatively large values of a, the original mode features a single-peak shape, while the new 
solutions are shaped as twin-peak modes, see Fig. [Hl^b). 

C. Bifurcations of the symmetric-symmetric (Sm-Sm) branches in the two-component system 

To understand the evolution of the Sm-Sm branches at finite a in the two-component system, we select a character- 
istic value, a — 0.6. The analysis was performed by fixing chemical potential A which is associated with component 4>, 
and scanning the configuration space of the system by letting vary the other chemical potential, ^, which is associated 
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Figure 8: (Color online) Saddle-node bifurcations which generate new pairs of modes with broken antisymmetry in the single- 
component model (word "single" in the plots stresses this fact) for a = 0.50 (a), a = 0.60 (b), a = 0.70 (c), and a = 0.80 (d). 
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Figure 9: (Color online) (a) A typical example of the pair of single-component modes with broken antisymmetry, which are 
generated by the saddle-node bifurcation shown in Fig. jS] This example is displayed for a — 0.60 and chemical potential 
|A| — 8.0 of the mode with the unbroken antisymmetry, (b) A typical example of the pair of symmetric single-component 
double-peak modes which are generated by the saddle-node bifurcation shown in Fig. 1101 This example is displayed for 
a — 0.90 and chemical potential |A| = 8.0 of the original mode with the unbroken antisymmetry. "Upper" and "lower" modes, 
generated by the saddle-node bifurcations, are defined according to the positions of the respective branches in Figs. ISland fTOl 



with component tp. At all values of A and /x, branches of solutions of the Sm-Sm types with unbroken and broken 
symmetries, have been found, the former ones always having lower energy ([S]), which suggests that the symmetry 
breaking may tend to destabilize the modes. 

The branch of the Sm-Sm solutions with the broken symmetry smoothly arises from the corresponding single- 
component asymmetric state, trapped in one of the nonlinear potential wells. In the course of its evolution following 
the variation of /i, this branch evolves towards the other well, where the other component of the "broken" Sm- 
Sm solution approaches the respective asymmetric state of the single-component model, while the mate component 
vanishes. 
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Figure 10: (Color online) The saddle-node bifurcation which generates a new pair of symmetric modes in the single-component 
model. At a = 0.60 (prior to the saddle-node bifurcation, panel (a)), only the "old" branch of the symmetric solutions is 
present. New symmetric branches, generated by the bifurcation, appear as a > 0.72. The new symmetric branches are plotted, 
along with the old one, for a = 0.80 (b), and for a — 0.90 (c). 



Similarly, the Sm-Sm branch with the unbroken symmetry arises from the corresponding symmetric solution of the 
single-component model, at some critical value of fx, and eventually merges into another symmetric state, in which, 
again, only a single component is present, but the other one. Thus, in both cases of the Sm-Sm modes with broken 
and unbroken symmetries, the evolution along the branch of the two-component solutions essentially amounts to 
the conversion of one single-component mode (resp., asymmetric or symmetric) into the mode represented by the 
originally missing component, while the opposite one vanishes. These bifurcation scenarios for the Sm-Sm branches 
arc illustrated by Fig. [TT] They persist in the interval of 0.16 < a < 0.72. 
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Figure 11: (Color online) The bifurcation scenarios, at a = 0.60, for the Sm-Sm modes with unbroken (a) and broken (b) 
symmetries, displayed by means of curves showing the evolution of energy ((6| with the variation of chemical potential ^, at the 
following fixed values of the chemical potential of the other component: [A] = 1.0,2.0,4.0,6.0 8.0. 



At a > 0.72, the system keeps the scenarios described above, but also gives rise to new ones, which involve the 
additional symmetric modes in the single-component model generated by the saddle-node bifurcation (see Figs. Iljb) 
andllOp. In particular, at a = 0.8 for fixed |A| = 4.0, there is a new branch of solutions with the unbroken symmetry, 
which smoothly arises from the lower symmetric state generated by the saddle-node bifurcation in the single-component 
model. As shown in Fig. I12f a). this two-component branch makes a loop in the norm space, {M, TV), and then merges 
into the upper symmetric state of the single-component model, which is generated by the above-mentioned saddle- 
node bifurcation. In this case, the second component of the Sm-Sm system vanishes at both ends of the branch. For 
larger absolute values of the chemical potential, say |A| = 8.0, a different behavior is observed. Instead of coming 
back to the upper state of the single-component model generated by the saddle-node bifurcation, the first component 
vanishes, while the other one evolves towards the corresponding upper-state solution of the single-component model, 
as shown in Fig. [T2T b) . 

As a increases further, |A| must be very large to support a link between the first and second components playing the 
role of the end states of the Sm-Sm branch. For a = 0.9, for example, both at |A| = 4.0 and |A| = 8.0, the Sm-Sm branch 
makes a loop in the norm space, emerging from the lower-state solution of the single-component model with the second 
component being absent. This branch evolves towards the upper-state symmetric solution of the single-component 
model with the second component vanishing again, as shown in Figs. [T^] (c,d). The further increase of a at relatively 
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Figure 12: (Color online) The same as Fig. 1111 but for |A| 
1.0 (the first, second, and third rows, respectively). 



= 4.0 and |A| = 8.0 (left and right columns), for a = 0.8, 0.9, and 



small values of |A| eliminates the looped branch; however, the branches of the two-component solutions with the 
broken and unbroken symmetries, which link the corresponding single- component modes (asymmetric or symmetric, 
respectively) still persist, as shown in Fig. [T2re). Nevertheless, at larger values of A the analysis reproduces the same 
pattern as observed before, cf. Figs. [El^d) and[T2Uf). 
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D. Bifurcations of the antisymmetric-antisymmetric (AS-AS) branches 

The structure of the configuration space of the AS-AS solutions in the model with finite a is complex. It can also 
be studied by fixing, say, a = 0.6 and one of the chemical potentials (A, which is associated with the ^-component), 
and scanning the space by varying /z, which is associated with component '0. In fact, the picture displayed below for 
a = 0.6 adequately represents the situation in the interval of 0.16 < a < 0.76. 

Similar to the case of the Sm-Sm modes considered above, the description of various bifurcation scenarios for the 
AS-AS branches, with the broken and unbroken antisymmetry alike, essentially reduces to identifying pairs of states 
of the single-component models which are linked by these branches. First, at relatively small fixed values of |A| (see 
Figs. [TST a.b)). we observe a branch of the AS-AS solutions with broken antisymmetry, which arises smoothly from 
the "unbroken" AS-AS family, at some critical value of /i. The asymmetry of the states rapidly increases along this 
branch. 

The branch extends up to a second critical value of fi, where it merges back into the "unbroken" AS-AS branch. 
Shortly afterwards, the unbroken branch itself merges into the corresponding solution of the single-component model. 
This scenario implies that, typically, one of the components, e.g., </>, gradually approaches the limit solution with the 
unbroken antisymmetry, while the other component remains strongly asymmetric almost until the end of the branch, 
where its asymmetry abruptly drops to zero. 

(a) (b) 



N-i 




(c) (d) 

AS-AS broken AS-AS broken 




Figure 13: (Color online) The picture of the bifurcation scenarios similar to that in Fig. 1111 but for AS-AS branches with 
a = 0.6 and A = 1.0 (a), 1.5 (b), 4.0 (c), and 8.0 (d). 

As |A| increases beyond a critical value, the above-mentioned broken-antisymmetry modes, generated by the saddle- 
node bifurcations in the single-component model, come into the play, making the evolution of the AS-AS branches 
quite involved. A new branch arises smoothly from the new broken- antisymmetry solution of the single-component 
model, which is trapped in one of the nonlinear potential wells. Fixing |A| = 1.5 and varying /i, we observe that this 
branch first makes a double loop in the norm space, (Af, N), with both components remaining trapped in the potential 
well, as shown in Figs. I13f c.d). Surprisingly, almost immediately another branch arises. The evolution along the new 
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branch amounts to a transition of one of the components from one nonUnear-potential weU to the other. One of the 
components on this branch looks hke an S-AS bound state, rather than a single AS state, as shown in Fig. [Ufb). 




Figure 14: (Color online) (a) and (b); Examples of AS- AS solutions belonging, respectively, to the lower (closed loop) and 
upper (transition) branches in Fig. I13f c). with a — 0.6, fi — 8.0 and |A| — 4.0. 

To consider in further detail how the evolution of the AS- AS branches depends on the value of a, we now fix chemical 
potential A and consider a range of values of a. In the case of not very large |A|, the AS-AS branch with the broken 
antisymmetry again originates from the "unbroken" AS-AS branch, as /i increases above some critical value. Since 
this branch is actually observed in the entire parameter space, including the limit case with the (5- functions (a — >■ 0), 
it may be naturally called the fundamental AS-AS branch with the broken antisymmetry. Along this branch, one of 
the components, say i^, rapidly grows with the increase of |^|, remaining localized in one nonlinear-potential well [as 
seen in Fig. llSf a)]. while the other component, 1p^ has two well-pronounced peaks near x = ±1. The evolution along 
this branch proceeds through the decrease of the maximum values of the second component, so that, at some point, 
the mode approaches the limit case of two AS states localized around opposite minima of the nonlinear potential, see 
Fig. llSf b). Finally, the system evolves towards the mirror image of the initial configuration, where the (j) component 
features maxima near x = ±1, while component ^ is trapped in the nonlinear potential well opposite to that where 
(j) was originally trapped, see Fig. llSf c). As fi approaches the second critical value, the branch merges back into the 
AS-AS state with the unbroken antisymmetry. Actually, this branch links two single-component AS solutions with 
the unbroken antisymmetry. 




''^ -3 -2 -1 1 2 3 4 "'^ -3 -2-1 1 2 3 4'^ -3-2-101234 

X " X 

Figure 15: (Color online) Examples of AS-AS solutions with the broken antisymmetry at three points of the fundamental 
branch with a = 0.8, |^| = 4.0 and |A| = 4.0. 



The evolution pattern becomes still more complex at larger values of |A|. In that case (for instance, at |A| = 4) , the 
building blocks of the AS-AS states are again single-component AS modes with the strongly broken antisymmetry, 
each localized around one minimum of the nonlinear potential. The evolution trajectory in the configuration space 
may then feature two loops. Along the first one, the system moves from the localization in one nonlinear potential well 
to the other, simultaneously swapping the two components (the one which was originally dominating and its initially 
vanishing counterpart). Continuing the evolution along the second loop, the system swaps the two components once 
again (not shown here in detail). 



E. Bifurcations of the symmetric-antisymmetric (S-AS) branches 

The evolution of branches of the S-AS type in the system with finite a also strongly differs from what was presented 
above for the system with the set of the (5-functions, corresponding to a — 0. As well as the solutions of the Sm-Sm 
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and AS- AS types, in the present case the evolution also amounts, essentially, to the switch between different states of 
the single-component model, linked by the branch of the S-AS type, with the unbroken or broken (anti)symmetries. 
To present the results, we again start with a = 0.6, fixing the chemical potentials (/x) which is associated with the 
symmetric component (</)), and scan the configuration space by varying the chemical potential (A) associated with the 
antisymmetric component, ip. 

Following the general pattern outlined above for other species of the two-component solutions, the "broken" S-AS 
branch arises from its counterpart with the unbroken (anti)symmetries, at some critical value of A. The energy of the 
emerging "broken" branch rapidly grows with |A|, as seen in Fig. 1161 
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Figure 16: (Color online) Branches of modes of the S-AS type: the evolution of the energy ([6]) in the norm space, {M,N), at 
a = 0.6 and = 0.2 (a), \n\ = 0.5 (b), \fi\ = 1.0 (c), and = 1.5 (d). 

The emergence and evolution of these S-AS branches is further illustrated by Fig. [TBI which displays the pattern 
of the modes at = 1.5,1.0,0.5,0.2. We observe that, if |^| remains relatively small, the increase of |A| leads to 
a bifurcation between the "broken" and "unbroken" S-AS modes. As |A| slightly decreases below this critical value, 
the norm of the AS component in the "unbroken" S-AS mode rapidly drops to zero, so that the branch almost 
immediately merges into the branch of symmetric solutions of the single-component model. As |A| increases above 
the critical value, the "broken" S-AS branch evolves towards the limit of a single-component "broken" AS solution 
which is strongly localized in one well of the nonlinear potential, while the "unbroken" S-AS branch approaches the 
"unbroken" AS mode of the single-component model. 

Another look at the evolution of the S-AS branches is suggested by the general circumstance, similar to that 
emphasized above, in the course of the analysis of the solution families of the Sm-Sm and AS- AS types, viz., that 
various branches of two-component solutions may be realized as links connecting the basic modes of the three types: 
the two-component ones with the unbroken (anti)symmetry, single-component symmetric/antisymmetric/asymmetric 
states, and the other above-mentioned species of single-component solutions, which are almost entirely localized in one 
well of the nonlinear potential. Below, we refer to these three basic types of the modes as A, B, and C, respectively. In 
particular, the picture of the evolution of the "broken" S-AS branches, as described above, may be interpreted so that 
the broken-(anti)symmetry branch evolves from state A towards B. For example, at small values of the "broken" 




M-norm 



20 



(a) 



(b) 



N 




0.8 



12 



2 



Figure 17: S-AS branches in the system with a = 0.60: the evolution of the asymmetry of the (originally) symmetric field 
(a) and antisymmetric field ^ (b) in the norm space, (M, A'^). 

S-AS states are weakly localized, and the increase of |A| drives the system towards state B without bifurcations. At 
> 0.15, the asymmetry of the AS component changes its sign, but there is still no bifurcation, with energy © 
monotonically depending on norms M and N. 

At > 0.37, there emerges the above-mentioned single-component mode (C), strongly localized in one well of 
the nonlinear potential, with a result that the AS state tends to be trapped near this mode, before it will be driven 
to mode B. This circumstance changes results of the bifurcation, as can be observed in terms of the asymmetries 
and energy in Figs. [T71 and [TW aV Another interpretation is that there is an additional state which appears at the 
bifurcation point, built as a combination of a symmetric mode in one component and a new AS mode with broken 
antisymmetry, see Fig [THl This means that three branches are observed in this region, with the evolving one attracted 
to state C and oscillating about it for a while, before plunging to the limit of the vanishing AS component (state B). 



Figure 18: (Color online) The evolution of modes of the S-AS type with the broken antisymmetry, approaching the effectively 
single-component mode (C, as defined in the text). In terms of Fig. [8l C is the upper AS state with the broken antisymmetry, 
generated by the saddle-node bifurcation and trapped in one nonlinear potential well. The parameters are a = 0.6, — 1.5, 
and |A| = 6.5 (a), 7.40 (b) and A = 7.43 (c), respectively. 

The picture of the evolution of the S-AS branches outlined above is valid in the interval of 0.1 < a < 0.67. Additional 
ramifications of the picture, which are not presented here, are revealed by the analysis at larger values of a. 



In this work, we have introduced the two-component one-dimensional model with the nonlinear (pseudo)potential 
represented by two strongly localized potential wells. The subject of the analysis is the evolution of various two- 
component modes of the three types: Sm-Sm (symmetric-symmetric), AS- AS (antisymmetric-antisymmetric), and 
S-AS (mixed states with symmetric and antisymmetric components), and the related SSB (spontaneous-symmetry- 
breaking) effects. In the limit of the nonlinearity modulation represented by the set of two (5-functions, the solution 
was obtained in the semi-analytical form. The consideration of this solution demonstrates new features of the two- 
component modes with spontaneously broken (anti)symmetries, which are qualitatively different from what was 
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Figure 19: (Color online) (a) The energy of the S-AS states with the broken (anti)symmetry at a = 0.6 and |A| = 
4.0,3.5,3.0,2.5,2.0,1.0,0.7,0.4, in the {M,N) plane, (b) Profiles of the two components of the new S-AS solution found 
in a vicinity of state C (defined in the text), at a = 0.6, and |/^| — 4.00, |A| = 0.791. 



previously reported in the single-component model. In particular, the spontaneous breaking of the antisymmetry is 
possible only in the two-component system, and, obviously, S-AS states, with the unbroken or broken mixed symmetry 
exist only in the two-component system. 

In the model based on the set of two nonlinear potential wells of a finite width, the evolution and bifurcation 
scenarios were found to be still more complex. First, new results for the case of finite potential wells were reported 
for the single-component model. These are pairs of modes with the broken antisymmetry, and pairs of twin-peak 
symmetric modes, generated by isolated saddle-node bifurcations, in either case. Then, numerous complex scenarios 
of the evolution of the Sm-Sm, AS-AS, and S-AS branches with the broken (anti)symmetries may be interpreted 
in terms of routes linking three species of basic modes: two-component states with the unbroken (anti)symmetries, 
"unbroken" or "broken" single-component states featuring density peaks in both nonlinear-potential wells, and, finally, 
single-component states effectively trapped in a single well. 

It remains to systematically analyze the stability of diverse two-component stationary modes reported in this work. 
Challenging problems may also be to extend the model to a set of several nonlinear spots, including a periodic lattice, 
and, on the other hand, to develop the two-dimensional (2D) version of the two-component system, following the 
recent analyses of the 2D extensions of the single-component model [l^, [2^ . 
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